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A suite of three evolution systems is presented in the 
framework of the 3+1 formalism. The first one is of second 
order in space derivatives and has the same causal structure 
of the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) system 
for a suitable choice of parameters. The second one is the 
standard first order version of the first one and has the same 
causal structure of the Bona-Masso system for a given pa- 
rameter choice. The third one is obtained from the second 
one by reducing the space of variables in such a way that the 
only modes that propagate with zero characteristic speed are 
the trivial ones. This last system has the same structure of 
the ones recently presented by Kidder, Scheel and Teukolski: 
the correspondence between both sets of parameters is explic- 
itly given. The fact that the suite started with a system in 
which all the dynamical variables behave as tensors (contrary 
to what happens with BSSN system) allows one to keep the 
same parametrization when passing from one system to the 
next in the suite. The direct relationship between each pa- 
rameter and a particular characteristic speed, which is quite 
evident in the second and the third systems, is a direct con- 
sequence of the manifest 3+1 covariance of the approach. 

I. INTRODUCTION 

In a recent article (which will be hereafter referred 
to as the KST paper), a wide class of hyperbolic first 
order formalisms has been studied with a view to Nu- 
merical Relativity applications. All these formalisms do 
use the well known 3 + 1 decomposition of spacetime 

ds 2 = -a 2 di 2 

+ 7^ {dx i + (3 l dt) {dx j + P j dt) i,j = 1, 2, 3 (1) 

which allows one to use Einstein's field equations as an 
evolution system for the metric ([!]), namely 

{d t - C^ij = -2 a K l3 (2) 
{d t - Cp)K iS = -V iaj (3) 
+ a [< 3) i?„- - 2K 2 + trK K l3 - Ry] 

(ADM system). The KST paper describes many ways of 
obtaining a first order version of the second order ADM 
system by introducing a number of arbitrary parameters. 
Different parameter choices correspond to different ways 
of using the energy and momentum constraints 

{3) R-tr{K 2 ) + {trK) 2 = 2t (4) 
V k K k i -d i {trK) = S i (5) 



to modify the structure of the resulting systems. The goal 
is to get well posed systems, so that existence, unicity and 
stability of the solutions can be ensured Q . In this sense, 
the KST paper can be understood as a generalization 
of previous works by other groups which have in 

common the fact that the only independent quantities 
are the metric coefficients and their first derivatives. 

A similar approach can be used to generalise the hyper- 
bolic first order formalisms proposed by the Palma group 
HQ, as suggested in Ref. M. The main difference with 
the formalisms studied in the KST paper is that in this 
case there are three supplementary variables, indepen- 
dent of the metric derivatives, whose evolution equations 
are obtained by reversing the order of space and time 
derivatives in the momentum constraint (JsJ) . This use 
of the momentum constraint to evolve three supplemen- 
tary variables can also be found in another context p|,p!o[ , 
where a second order system of the form (||) is considered 
(BSSN system). 

To summarize, we have different formalisms in which 
the constraints (Q, ||) are used to modify the structure of 
the ADM evolution system. The momentum constraint 
is used in some cases to evolve supplementary quantities, 
either keeping the system to be of second order (|jl^] or 
getting first order versions But the same momen- 

tum constraint is used instead in other cases to modify 
the first order versions without introducing extra quan- 
tities g-§. 

The purpose of the present article is to point out the 
strong relationship among all these formalisms from the 
theoretical point of view. With this aim, we propose a 
new covariant framework in which the momentum con- 
straint is used to evolve a vector quantity, so that the 
properties of the resulting systems are independent of 
the choice of space coordinates. 

The paper is organised as follows: In Section II, we in- 
troduce the vector quantity Zi, which will replace the 
corresponding non-vector supplementary quantities in 
0-jni- The original ADM system is then generalised by 
making use of the vector Zi and the energy constraint (Q) . 
In Section III, the resulting second order system (which 
will be hereafter referred to as system A) is compared 
with the BSSN system [gjilj. 

In section IV, we obtain a first order version of sys- 
tem A (which will be hereafter referred to as system B), 
which is compared with the formalisms presented in ^|,[7| . 
Finally, in section V, system B is further modified in or- 
der to be compared with the formalisms discussed in the 
KST paper. 
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II. INTEGRATING THE MOMENTUM 
CONSTRAINT 

Let us define the vector Z % as follows: 

(ft - C p )Z t = a [Vj (Ki j - Si j trK) - St] (6) 

so that it is clear that Zi vanishes for physical solutions, 
which can be defined as the ones verifying both the en- 
ergy and momentum constraints ([|, ||). The "zero" vec- 
tor Zi can then be considered as a good covariant measure 
of the deviation of any solution (a, jij, Kij) of the ADM 
system from the submanifold of physical solutions, as far 
as the momentum constraint is concerned. This means 
that the cumulative effect of the differential constraint 
(||) can be allowed for by looking at the simple algebraic 
equation 

Zi = (7) 

which is to be regarded as the integral version of (|J). 

But these Zi, besides being useful auxiliary quanti- 
ties, can be also considered as supplementary dynami- 
cal variables. This is the kind of approach proposed in 
Rcf. jy], where a "A-system" was proposed containing 
as much as 40 supplementary quantities to be evolved 
by using the constraints. We will limit ourselves here to 
only three quantities Zi (corresponding to Xi in Jll|]), the 
ones related with the momentum constraint. We will use 
these quantities Zi, together with the energy constraint 
to implement covariant quasilinear modifications of 
the principal part of the ADM system, without affecting 
the submanifold of physical solutions (the ones verifying 
(H @))- We will call "system A" the most general second 
order evolution system obtained in this way, namely 

(d t - C P )K ZJ = -ViOtj (8) 
+ a [WRij - 2Kl + trK K %] - 

-^-[R-tr{K^) + tr^K)-2r] llJ 
+ -fWiZj + VjZi) - — (V • Z) 7ij 

where v and n are arbitrary parameters. 

The modified evolution system (||,||) (system A) is of 
a mixed type: first order in (Zi,Kij) but second order 
in (a,jij). In this sense, it is very similar to the BSSN 
system In order to compare these systems, let us 

state some definitions: 

a) Two evolution systems will be said to be equivalent 
if one can be obtained from the other by a change of 
variables. It is clear that two equivalent systems have 
the same space of solutions. 

b) Two evolution systems will be said to be quasiequiv- 
alent when their principal parts are equivalent. It follows 
that two quasiequivalent systems have the same causal 
structure. 



For instance, if we introduce "damping terms" to mod- 
ify eq. (^) (as proposed in Jll[] with a view to Numerical 
Relativity applications), we are not affecting the princi- 
pal part (damping terms are algebraic in Zi), so that the 
resulting "damped" system would be quasiequivalent to 
system A. Also, if we introduce additional terms of the 
form ZiZj or Z 2 ~/ij in (^), the resulting modified system 
would be quasiequivalent to system A as well. 

On the contrary, different values of the parameters n, 
v, n in (|^) do correspond to different forms of the prin- 
cipal part, so that they would lead to different evolu- 
tion systems which are not quasiequivalent, even if they 
have (as they do) the same submanifold of physical so- 
lutions. In this context, it is natural to use system A as 
a framework to compare different formalisms that have 
been proposed with a view to Numerical Relativity ap- 
plications, where the causal structure of the evolution 
system is important to ensure the well posedness of the 
problem, namely the existence, unicity and stability of 
the solutions. 

We can also take advantage of the covariance of system 
A under coordinate transformations of the generic form: 

t' = f(t) (9) 
x 1 ' = g{x\t) 

(3+1 covariance) so that we can choose j3 l = (normal 
coordinates) in what follows without any loss of general- 
ity. 

III. COMPARING WITH THE BSSN SYSTEM 

Let us consider here the system introduced in by Shi- 
bata and Nakamura ||] and Baumgarte and Shapiro 
(BSSN system). The metric coefficients 7^ are expressed 
in terms of a conformal metric: 

%=e- 40 7y (10) 
with unit determinant, so that 

e 4 ^ 7 1 /3 = [de t( 7y . ) ]l/3 (11) 

The second fundamental form Kij is then decomposed 
into its trace and trace- free components, namely 

K = 7 « K^ (12) 

A 3 =e-^ {Kij-\ Km) (13) 

so that the principal part of their evolution equations can 
be written as p0| : 

d t K = - 7 « V, aj + ... (14) 
d t Aij = e~ 4 * (-V< aj + a R lJ ) TF + ... (15) 

where the superscript TF denotes the trace-free part of 
a tensor. 
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The principal part of the Ricci tensor in ( |15| ) is decom- 
posed as follows [ [Tof : 



-2d 2 



(16) 



where the "conformal connection" quantities T l can be 
defined as: 



(17) 



Up to now we have only performed a transformation of 
variables of the original ADM system. But notice that 
the system (14, [lq) contains seven evolution equations, 



compared with the six equations in ([|) . The extra quan- 
tity is the trace of Ay, which is supposed to vanish by 
construction. This opens the way to two different ap- 
proaches: 

a) to consider tr(A) as an independent variable, which 
can take non-zero values during numerical evolution. The 
equation 



tr(A) = 



(18) 



is to be considered as an additional constraint which can 
be used to monitor the accuracy of the simulation (free 
evolution approach). 

b) to enforce the vanishing of tr(A), by removing sys- 
tematically any contribution that could arise from trun- 
cation errors in the discretization of equation (15), so 
that tr(A) is not considered as a new dynamical variable 
(constrained evolution approach) 

In what follows, we will adopt the second approach, so 
that the system (J14|, |l^) is then quasiequivalent to the 
ADM system. 

The BSSN system does consider the P as independent 
dynamical variables, so that the full set of independent 
BSSN variables is given by 



{a, 0, 7y, K, Ay, f 1 } 



(19) 



where we must remember that both the trace of Ay and 
the determinant of "fy do have fixed, non-dynamical, val- 
ues. An evolution equation for the T l could be derived 
fiof by permuting a time derivative with the space deriva- 
tive in @ 



d t T l = -2 a A 1 ; 



(20) 



But the BSSN system uses instead the momentum con- 
straint (|]) to transform eq. (^0|) into 



(21) 



Now we are in position to compare the BSSN system 
with the one introduced in the previous section (system 
A), with a set of variables given by 



= {<*, 7ij, Kij, Zi} 



(22) 



The conformal decomposition (10, 53], O, |T^) can be eas- 
ily performed on system A as well. The key point is then 
the passage from ( p0| ) to the BSSN evolution equation 
(pl|). In the framework of system A, the momentum con- 
straint is just the right-hand-side of the evolution equa- 
tion for Zi. It follows that the use of the momentum 
constraint to transform (the principal part of) the right- 
hand-side of (^) into (|2l|) does correspond to the follow- 
ing transformation of T 1 : 



~kj 

-Jtk 7 j 



2Zi 



(23) 



which is obviously consistent with the previous expres- 
sion ([it]) if we remember that Zi = for the physical 
solutions. 

Now we can compare easily (|l^) with the trace-free 
part of (H): it follows that they coincide if and only if 
[i = 2. The same comparison can be performed again 
between ( |l4] ) and the trace part of (||): it follows that 
they coincide if and only if v = n = | . 

To summarize, we have shown then that the BSSN 
system is quasiequivalent to system A for the parameter 
choice 



M — 2, v 



4/3 



(24) 



provided that the "conformal connection" quantities Y % 
in the BSSN system are related with the vector Zi by 
equation (p3|). 



IV. A FIRST ORDER EVOLUTION SYSTEM 

A first order version of the system A can be obtained 
in the standard way by considering the first space deriva- 
tives 



A k = d k (lna) 

D klJ = - <9fc7y 



(25) 
(26) 



as independent dynamical quantities whose evolution 
equations are given by 



d t A k + d k (a Q) = 
d t D kij + d k (a Ky) = 

where the quantity Q, defined as 

dt In a — —a Q 



(27) 
(28) 



(29) 



can be related with trK in order to fix the time coordi- 
nate gauge, namely 



Q = f trK 



(30) 



where / is an arbitrary function of the lapse a. 

This allows to transform system A into a first order 
evolution system for the set of variables 
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{a, jij, Kij, Ai, D ktj , Zi] 



(31) 



with a non-trivial principal part given by M, B7j, |2q) and 



dtK i3 +d k (a A%) 
where we have noted 



\ k - 



D k tJ + ^iiAj +Dj-2 D rj r - f i Zj) 

+ -S k j(Ai + A - 2 D ri r - fiZi) 
-\{D k -W k r)ll] + V -Z k 7ij 



(32) 



(33) 



and 



D k = Y 3 D 



kij 



(34) 



This first order system will be called system B in what 
follows. 

System B is very similar to the Bona-Masso system 
(we will follow here the notation of Ref. Q] to avoid 
confusion) , which is expressed in terms of the set of vari- 
ables 



u = {a, 7y, K^, Ai, Dkij, Vi} 
where the quantities Vi can be defined as 

Vi = Di - Dji j 



(35) 



(36) 



An evolution equation for the Vi could be derived by 
combining (36) with the evolution equations ( |28| ) for the 

Dkij 



d t Vi = adj {Ki> -trK 8f) + ... 



(37) 



but the Bona-Masso system uses instead the momentum 
constraint to transform ( |37| ) into an evolution equation 
with a trivial principal part, namely 



dtVi 



(38) 



The remaining evolution equations in the Bona-Masso 
system are identical to those in system B, with the re- 
placement of (B3h by 



(39) 



\ k ij =D k ij + ~8 k i (Aj + 2Vj-D j ) 
+ ls k j (A i + 2Vi-D i )-^V k lij 



In order to compare system B with the Bona-Masso sys- 
tem, the key point is again the passage from ( p7|) to the 
Bona-Masso evolution equation ( |38|) . The use of the mo- 
mentum constraint to transform (the principal part of) 
the right-hand-side of ( |37| ) into (pq)) does correspond to 
the following transformation of Vc 



V, = A - »,. ' - Zi 



(40) 



which is again consistent with (|3 
for physical solutions. 



because Zi vanishes 



Now we can substitute ( 40 ) into ( |39| ) and compare with 
([$||): it follows that they coincide if and only if 



(41) 



To summarize, we have shown that the Bona-Masso sys- 
tem is quasiequivalent to system B for the parameter 
choice (|4l|), provided that the supplementary quantities 
Vi in the Bona-Masso system are related with the vector 
Zi by equation (p0[). As a consequence of this fact, sys- 
tem B inherits (for the parameter choice (|4l"|)) the causal 
structure of the Bona-Masso system, which is known to 
be strongly hyperbolic (real characteristic speeds and a 
complete set of eigenfields 0j) if and only if / > 0. 



V. REDUCING VARIABLES SPACE 

System B has two sets of non-trivial standing eigen- 
modes (the ones with zero characteristic speed in normal 
coordinates). The first one is related with gauge evolu- 
tion, namely 



Ai- f Di 



(42) 



while the second one is the combination ( |4C| ) , related with 
the momentum constraint. We will see how we can use 
this fact to reduce the set of independent dynamical vari- 
ables ([H]) of the system B to the smaller one used in the 
KST paper, namely 



u = {yiji Kij, dkij} 
where, in spite of the fact that 



dkij — 2 D k ij 



(43) 



(44) 



we will keep using here the KST original notation to avoid 
confusion. 

First of all, one can explicitly integrate the relationship 
( p9 , 30) between the lapse and the volume element in the 
case / = 2 a = constant, namely 



d t (a 7 - ff ) = 



(45) 



so that the value of a can be defined in terms of 7 for 
every initial condition (densitized lapse) . The same thing 
can be done with Ai and Di , so that 



Ai = 2 a Di 



(46) 



The list of dynamical variables is then reduced from ( |3l| ) 
to 



= {lij, K i: j, D ki j, Zi\ 



(47) 



where the gauge-related standing modes ( f42| ) have just 
disappeared from the dynamical system. 
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In doing this, we have restricted ourselves to the 
/ = constant case, in spite of the fact that this case 
is known fO] to be prone to numerical gauge instabilities 
if / 7^ 1. We do that with a view to mimic the KST 
system, where a is defined to be a constant and, as a 
consequence of it, only examples with a = 1/2 (/ = 1) 
are provided. In constrast, system B would allow more 
general gauge choices, like the "1+log" one (/ = 2/ a) 
which is frequently used in Numerical Relativity. 

Disposing of the supplementary variables Zi in ( |47| ) is 
not so easy. We will take as a guide the KST evolution 
equations for 7^ and K^, which are very similar to the 
corresponding ones in system B: we need just to replace 



9 \ k ■ — rl k ■ ■ 

% 3 l J 



^i[(l + 2a) d,--2dr/ 



(48) 



-5 k 3 [{\ + 2a)d t -2 d ri r ] - - (d k - d rk r ) ll3 



where we have taken £ = — 1, 7 — —n/2 in the original 
KST expression to reproduce (|33| ) more closely. 

The main difference between system B and the KST 
one comes instead from the evolution equations for dkij , 
namely 

d t dkij = -2 a d k K i:j + x a 7^ (d r K r k - d k trK) (49) 
+ | a[ lkl {d r K r - djtrK) + lkj (d r K r i - dJrK)] + ... 

where n, x are arbitrary parameters. These equations are 
obtained by using the momentum constraint (|5|) to mod- 
ify the standard equations ( 28 ) which were used in system 
B. The key point again is to realize that the use of the 
momentum constraint to transform (the principal part 
of) the right-hand-side of ( |28| ) into that of (E9h amounts, 
allowing for (||), to modify the relationship (|44|) between 
D k ij and d k ij in the following way: 



2 Dkij = d k ij - V 7fe(; Zj) - x Zk-fij 



(50) 



which is consistent again with (|26|) because Zi vanishes 
for physical solutions. A straightforward substitution of 
( |46| ) and (|5(]) into ([33]) shows that the vector Zi disap- 
pears from the principal part if we just relate the free 
parameters in system B with those of the KST system as 
follows: 



/ = 2 a 

v = X - n(x - 1)12) 

n = V - X/ 2 - cr(v + 3 X) 



(51) 



so that the Einstein-Christoffel system jj| (corresponding 
to a = 1/2, 7} — 4, 7 = x = in the KST paper) is 
recovered when 



/ = 1, /i = 2, ^ = rt = 



(52) 



To summarize, we have shown here that the KST sys- 
tems are quasiequivalent (in the case £ = — 1) to a first 



order system which can be obtained by reducing dynam- 
ical variables in system B to the minimal set ([43|). The 
reduced system will be called system B' in what follows. 
The four parameters of system B' can be expressed in 
terms of those of the KST system by (|5l|), where we 
must remember that n = —2 7. 



VI. CONCLUDING REMARKS 

The evolution systems A, B and B' presented in this 
paper are closely related in a transparent way, up to 
the point that we have been able to keep the same set 
of parameters (/, fj,, v, n) in the three cases. The fact 
that, for suitable parameter choices, they are quasiequiv- 
alent respectively to the BSSN §,0, Bona-Masso §0] 
and Einstein-Christoffel || systems points out the strong 
relationship between these systems, proposed indepen- 
dently in different contexts by different authors, all of 
them extremely useful in Numerical Relativity applica- 
tions. 

There are at least two specific points in which the suite 
A, B, B' can improve our understanding of the previous 
systems. The first one is that the supplementary quanti- 
ties Zi in systems A and B are three-dimensional vectors, 
contrary to what happens with their counterparts T\ Vi 
respectively, which do not have tensor behaviour. Then, 
systems A and B can be used to build up covariant coun- 
terparts of the second order BSSN [p|,[io| and first order 
Bona-Masso p,0 systems, respectively. 

The second point can be seen when analysing the list 
of non-trivial characteristic speeds in the KST systems, 
namely 



{±1, ±ci, ±c 2 , ±c 3 } 
where in our case (£ = —1) we have j 



2 -2a 



(53) 



(54) 



C2 



C3 



^(2n-2ncr-x-6xv) 
1 + 2 7 - 7 77 + x + 2 7 x 



which can be translated, allowing for (J5lJ) into 

ci 2 = .f, c 2 2 =/i/2, c 3 2 = l + ^-n (55) 

where we have used n = —2j. 

Expression ( |55| ) allows to decouple the parameter set in 
a way such that every parameter is related with just one 
characteristic speed, making the physical interpretation 
more transparent: the gauge parameter f is associated 
only with c\ (gauge speed), the parameter u (which ap- 
pears in the terms affecting to the longitudinal part of 
Kij) is associated only with C2, whereas the two param- 
eters v and n appearing in the transverse trace terms 
are related with C3 (notice that there is some degree of 
redundancy between v an n). 
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It is not surprising then that in all the particular sets of 
parameters corresponding to previous works PHIo| (|24| , 
[II]. |5^) one has 

jit = 2, v = n (56) 

which amounts to the "physical speed" requirement for 
the degrees of freedom not related with the gauge (C2 2 = 
C3 2 = 1). We are aware that in the case of the BSSN 
system we can not speak properly about characteristic 
speeds, but the close relationship between systems A and 
B can bring some light even in this case (see for instance 
Ref. [^3) for a similar approach). 

The question arises of wether systems A, B, B' will 
perform any better than their well known counterparts 
(BSSN, Bona-Masso, KST) in numerical simulations. We 
have done preliminary tests, comparing system B against 
the Bona-Masso one in black-hole spacetimes. Our re- 
sults indicate that the Bona-Masso code performs bet- 
ter, because of the fact that the additional quantities Vi 
evolve with an ordinary differential equation j3q) , avoid- 
ing the truncation errors inherent to the discretization of 
the partial derivatives that appear instead in the evolu- 
tion equation (^|) for Zi. The fact that Zi is a true tridi- 
mensional vector does not help very much in that con- 
text. Something similar happens when comparing sys- 
tem A with BSSN: the evolution equation ( pl| ) for the 
T l is still simpler than (^|). In our opinion, the BSSN, 
Bona-Masso and KST systems are the state of the art 
(each one in its class) and will be difficult to beat by any 
other quasiequivalent system. One would need to pro- 
vide something more different, modifying even the causal 
structure of the system (which is given by the principal 
part), in order to improve in a significant way the perfor- 
mance (stability and accuracy) of the current numerical 
simulations. 
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